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The existence of a sizable, O (l0 _10 -10 _9 G) , cosmological magnetic field in the early Universe has 
been postulated as a necessary step in certain formation scenarios for the large scale 0([iG) magnetic 
fields found in galaxies and galaxy clusters. If this field exists then it may induce significant mixing 
between photons and axion-like particles (ALPs) in the early Universe. The resonant conversion of 
photons into ALPs in a primordial magnetic field has been studied elsewhere by Mirizzi, Redondo 
and Sigl (2009). Here we consider the non-resonant mixing between photons and scalar ALPs 
with masses much less than the plasma frequency along the path, with specific reference to the 
chameleon scalar field model. The mixing would alter the intensity and polarization state of the 
cosmic microwave background (CMB) radiation. We find that the average modification to the CMB 
polarization modes is negligible. However the average modification to the CMB intensity spectrum 
is more significant and we compare this to high precision measurements of the CMB monopole 
made by the far infrared absolute spectrophotometer (FIRAS) on board the COBE satellite. The 
resulting 95% confidence limit on the scalar-photon conversion probability in the primordial field 
(at 100 GHz) is P 7++( £ < 2.6 x 10~ 2 . This corresponds to a degenerate constraint on the photon- 
scalar coupling strength, g B g, and the magnitude of the primordial magnetic field. Taking the 
upper bound on the strength of the primordial magnetic field derived from the CMB power spectra, 
B\ < 5.0 x 10 _9 G, this would imply an upper bound on the photon-scalar coupling strength in the 
range g s < 7.14 x 10" 13 GeV _1 to g e « < 9.20 x 10 _14 GeV _1 , depending on the power spectrum of 
the primordial magnetic field. 

PACS numbers: 



I. INTRODUCTION 

The origin of the observed large-scale magnetic fields of 
order ^G found in nearly all galaxies and galaxy clusters 
is still largely unknown. It is generally believed that the 
galactic magnetic field develops by some form of ampli- 
fication from a pre-galactic cosmological magnetic field. 
The two popular formation scenarios that are considered 
are either some exponential dynamo mechanism which 
amplifies a very small seed field of order 10~ 30 G as the 
galaxy evolves, or the adiabatic collapse of a larger ex- 
isting cosmological field of order (10 -10 -10 -9 )G. There 
are a number of pros and cons to both scenarios. See 
[l| for reviews on the subject. Theories explaining the 
origin of this primordial magnetic (PMF) field are still 
highly speculative. It has been suggested that a large- 
scale magnetic field could be produced during inflation 
if the conformal invariance of the electromagnetic field is 
broken; see for example 2] or more recently Q. 

To date, there is no astrophysical evidence for the ex- 
istence of a large-scale cosmological magnetic field, and 
only upper bounds on its magnitude have been derived. 
So far the strongest constraints have come from measure- 
ments of the cosmic microwave background (CMB) and 
big-bang nucleosynthesis [1, 4-7]. The CMB bounds are 
derived by considering the effects of Faraday rotation in- 
duced by the PMF on the CMB power spectra. In [5] an 
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upper limit in the range 6 x 10 -8 to 2 x 10 _6 G was derived 
for the mean-field amplitude of the PMF at a comoving 
length scale of 1 Mpc by comparison to the WMAP 5- 
year data. More recent work [y, 0] analysing the WMAP 
data in combination with other CMB experiments such as 
ACBAR, CBI and QUAD have placed tighter constraints 
on the amplitude of the PMF with an upper bound on 
the mean-field amplitude at 1 Mpc of ~ 5 x 10 _9 G. 

The existence of a primordial magnetic field would in- 
duce mixing between CMB photons and axion-like par- 
ticles (ALPs). ALPs refer collectively to any very light 
scalar or pseudo-scalar with a linear coupling to F ilv F p,v 
or e-tivpoF^" F pa respectively. In this paper we consider 
non-resonant mixing of scalar ALPs and CMB photons, 
with specific reference to the chameleon scalar field model 
[8l4ic|. Standard ALPs have constant mass and photon- 
scalar coupling everywhere, while the chameleon model 
has a density dependent mass. In sparse environments, 
such as the primordial plasma, the chameleon acts as 
a very light scalar field and would be indistinguishable 
from a standard scalar ALP. However in dense envi- 
ronments the chameleon is very heavy and evades the 
standard ALP constraints The current best con- 

straints on the coupling strength between the chameleon 
and electromagnetic fields are: g e s < 9.1 x 10 _10 GeV _1 
[l(j and g cS < (0.72 - 22) x Hr 9 GeV _1 (T^. Other 
chameleon-like theories exist such as the Olive-Pospelov 
model [H, [l3| which have a density-dependent coupling 
strength, but we do not discuss these further here. 

The resonant mixing of photons with scalar and 
pseudo-scalar ALPs in a primordial magnetic field has 
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been analysed by Mirizzi, Redondo and Sigl [l4|. They 
find a constraint on the combined magnetic field strength, 
B, and ALP-photon coupling strength, g: g{B 2 ) 1 / 2 < 
lCT 13 - 10- 11 GeV" 1 nG, for ALP masses between 
10~ 14 eV and 10 _4 eV. We believe that these results will 
not necessarily be applicable to the chameleon because 
its mass evolves as the density of the Universe decreases. 

In the following analysis we assume a stochastic pri- 
mordial magnetic field with a power-law power spectrum, 
similar to the treatment in Q. We assume fluctuations 
in the magnetic field are damped on small scales due 
to Alfven wave dissipation [l5j], and subdivide the mag- 
netic field into multiple domains of length L of a compa- 
rable size to just above the Alfven wave damping scale 
kp . The magnetic field in each domain is assumed to 
be approximately constant, and correlated to the other 
domains according to the magnetic power spectrum. A 
similar method was applied to correlations in quasar po- 
larization spectra by Agarwal, Kamal and Jain [r| . 

The degree of conversion between chameleons and 
photons in a magnetic field is inversely proportional 
to the electron density in the plasma [l^- Hence the 
dominant contribution to photon-scalar mixing will take 
place in the region after recombination when the ioniza- 
tion fraction drops significantly and before reionization. 
This greatly simplifies the mixing equations because we 
do not need to evolve the photon-scalar mixing equa- 
tions through the last scattering surface, nor include the 
density inhomogeneities present after reionization. We 
model a scenario in which the primary CMB is formed 
at the last scattering surface and then evolves through a 
primordial magnetic field extending from recombination 
(z ~ 1100) to the epoch of reionization (z ~ 20). 

This paper is organized as follows: in section |TT] the 
chameleon model is introduced and the calculations de- 
scribing photon-scalar mixing in a magnetic field, living 
in a Friedman-Robertson- Walker (FRW) spacetime, are 
presented. The power spectrum of the primordial mag- 
netic field is discussed in more detail in section Mil In 
section IIVI we analyse the evolution of the photon and 
chameleon states as they propagate through the multiple 
magnetic domains, and predict the average modification 
to the CMB intensity and polarization. In section [Vj our 
predictions are compared to precision measurements of 
the CMB monopole made by the far infrared absolute 
spectrophotometer (FIRAS) on board the cosmic back- 
ground explorer (COBE) satellite. We present a sum- 
mary of the work and our conclusions in section IVD The 
appendix contains details of the equations governing the 
evolution of the CMB Stokes parameters. 



II. CHAMELEON-PHOTON MIXING IN AN 
EXPANDING UNIVERSE 

The chameleon scalar field has been suggested as a can- 
didate for the dark energy quintessence field Q. Stan- 
dard quintessence scalar fields have coupling strengths 



to matter which are unnaturally fine-tuned to very low 
values so as to be compatible with fifth-force experi- 
ments. By contrast the chameleon model is constructed 
so that the chameleon can have a gravitational strength 
(or greater) coupling to normal matter while at the same 
time evading fifth-force constraints. This is achieved 
by introducing a density-dependent term in the effec- 
tive potential of the chameleon, which causes a change 
to the minimum of the potential depending on the den- 
sity of the surrounding matter. In sparse environments 
the chameleon behaves as an effectively massless scalar 
field while in a laboratory on Earth it is much heavier 
and evades detection. 

In addition to the matter coupling, the chameleon can 
have a non-zero coupling, j^—F^F^, to the electro- 
magnetic (EM) field. The mass parameter, M e g, de- 
scribing the strength of the photon coupling is best con- 
strained by observations of radiation passing through as- 
trophysical magnetic fields. The two best current con- 
straints in terms of the coupling strength g c g = 1/M e g 
were presented in section [IJ In terms of M e g, the lower 
limit from considering constraints on the production of 
starlight polarization in the galactic magnetic field is 
M cff > 1.1 x 10 9 GeV Measurements of the Sunyaev- 
Zel'dovich effect in galaxy clusters places a lower bound 
on M off in the range 4.5 x 10 7 GeV to 1.4 x 10 9 GeV, de- 
pending on the model assumed for the cluster magnetic 
field [f|. 

The action describing the chameleon model is that of 
a generalized scalar-tensor theory: 

S = J d^x^-g {^M 2 P {R - ^ v d^d v <t> 
-V(cj>) - ^B F { < t>/M)F ltv F'A 

+Smattcr (V W ,ff$) > 

where V(4>) is the self-interaction potential of the scalar 
field 4>, and <S ma ttcr is the matter action, excluding the ki- 
netic term of electromagnetism, which contains the cou- 
pling of the matter fields ip^ to the conformal met- 
ric g^l = Bf(<f>/M)g ftv . The B t ((f>/M) determine the 
coupling of the scalar field to different matter species, 
and Bp((f>/M) determines the photon-scalar coupling. 
For simplicity a universal matter coupling, Bi(<f>/M) = 
B m (4>/M), is assumed. We take the Universe to be de- 
scribed by a spatially-flat Friedman-Robertson- Walker 
(FRW) metric: ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) 
where a(t) is the time-dependent scale factor describing 
the expansion, normalized to ao = 1 today. The coor- 
dinates (x, y, z) are comoving, related to the physical 
coordinates by and t is the proper time. 

The equation of motion for the <f> field is found by vary- 
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ing S with respect to 



n<t> 



3H{t)d t 



V 2 < 



B i(4>) B 'i(<f>)Tmi 



where H(t) = a/a is the Hubble expansion rate. T$ 
is the trace of the stress-energy tensor in the conformal 

(i) 

frame described by the metric g)il, for which particle 
masses are constant and independent of <\>. This is re- 
lated to the stress-energy tensor in the physical frame 

by T m = BfT®. In general, T m = -{p m - 3P m ) corre- 
sponding to the physical, measured density and pressure. 
If we assume the chameleon does not couple to the en- 
ergy density of dark matter or any other exotic particles, 
then T m — —p h , the baryonic matter density. Hence 



□0 = K' ff (^; F 2 ,Pb), 



where, 



V eS {cj>\ F 2 , Pb ) 



V(<t>) + -F 2 B F {<jy) + p b B m {cj>). 



The form of the self-interaction potential V(<f>) de- 
termines whether a general scalar-tensor theory is 
chameleon- like or not. For a chameleon field we require 
the potential to be of runaway form. A typical choice of 
potential can be described by 



V(<f>) w A + 



A™+ 4 
n(j} n 



(1) 



where for reasons of naturalness A ~ O (Aq). If the 
chameleon is to be a suitable candidate for dark energy, 
we require A = (2.4 ± 0.3) x 10 _3 eV. The chameleon 
mass is defined by 



vZ (0 mi n; f 2 , Pb ) 



where F^ v is the background value of the electromagnetic 
field tensor. 

The photon equation of motion comes from varying S 
with respect to SA^: 

V v (B F (<f>)F> w ) = 0. 

Any electromagnetic components in L ma tter would intro- 
duce electromagnetic currents on the right-hand side of 
the equation. We neglect these at this stage in the calcu- 
lation. The primordial plasma is a good conductor and 
so any currents will be small [17]. To a good approxima- 
tion the induced currents arising from photon propaga- 
tion will be described by the plasma frequency which is 
included later in this calculation. 

In standard chameleon theories we assume <f>/M -C 1 
and approximate Bp ~ 1 + 4>/M c s and B 2 ^ s» 1 + 2<p/M, 
where this defines M c g and M. We expect these two 
coupling strengths to be of similar magnitude but do not 



require it. Under this approximation, and splitting the 
electromagnetic field into the background EM field F^ u 
and the photon field f^, the above equations of motion 
become 



Uip ~ m%<p + -— (U„F^ + rP^) 



4M t „ 



where tp = cj> — <f> is the perturbation in the scalar field 
about its background value, and we neglect terms that 
are O (f 2 ) and 0(<Pf» v /M eS ), 

In an inertial, locally Minkowskian frame, with metric 
9ij,u = diag (— 1, 1, 1, 1), the electromagnetic field tensor 
has the form 

E x E y E z 
piiv _ | -F x B z -By 
-E y —B z B x 
-E x By —B x 

The coordinate transformation from the locally 
Minkowski metric, g^ — diag (— 1, 1, 1, 1), into the 
FRW metric, = diag (—1, a 2 , a 2 , a 2 ), is 



with 



K = = diag (1,1/a, l/o,l/o). 

Thus in an FRW expanding Universe the EM field tensor 
is given by 



ptiv 



E x /a Ey/a E z /a 

-E x /a B z /a 2 -B y /a 2 

-Ey/a -B z /a 2 B x /a 2 

-E z /a By/a 2 -B x /a 2 



We take the background to be a large-scale magnetic 
field B, and assume the contributions from electric fields 
in the plasma are negligible. The photon field can be 
described by the quantum polarization states a M where 
/ M „ = V^a^ — V„a M . Taking the Lorentz gauge condi- 
tion V M a M = d^a^ 1 + 3Ha° — 0, we find the equations of 
motion become 

-d 2 tV - 3H(t)dt<p + -^VV = m%<P + B ' (V X a) 



M 



off 



and 



„n „ n n 1 n V<j9 X B 

-d?a.-5Hd t a + 2qH 2 a-4H 2 a.+ — V 2 a = — , 

a 1 a l M c R 

where we define the deceleration parameter, q(t) = 
—da/a 2 . 

In this analysis we have assumed that the background 
values of the chameleon and photon fields are slowly vary- 
ing over length and time scales of 0(1 /lS) where to is the 
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proper frequency of the electromagnetic radiation being 
considered. We further assume that the frequency of the 
radiation satisfies uj >• H and that q < 0(1). The equa- 
tions describing mixing between the photons and scalar 
field are then 



1 

-a+ — V 2 a 

a 1 



Vif x B 

a 2 M, 



cir 



B ■ (V x a) 



off 



where we have included the plasma frequency, 
4ira em n c /m c , as an effective photon mass. Whenever 
electromagnetic radiation propagates through a plasma 
with electron number density n c it displaces the elec- 
trons slightly and induces oscillations at their natural 
frequency uj p \. This interaction of the EM wave with 
the electron density hinders its progress and acts as an 
effective mass for the photons. 

Taking the radiation to be propagating in the z direc- 
tion of an orthonormal Cartesian basis (x, y, z), with a = 
(lx,1y,0) T , the equations of motion for the chameleon 
and photon can be written in matrix form: 



d 2 



V 



"pi 


A/off 







pi 



Byd Z 

a 2 M ctt 
Bxd z 



m 



\7y) 



= 0. 



Note how the chameleon scalar field only mixes with the 
component of photon polarization aligned perpendicular 
to the transverse magnetic field. The magnitude of the 
magnetic field aligned parallel to the photon path (B z ) 
plays no part in the mixing. 

Following a similar procedure to that in [l2|, [l8| we 
assume the fields vary slowly over time compared to the 
frequency of the radiation and that the refractive index 
is close to unity, which requires m'i/2oj 2 1 ui 2 i /2uj 2 and 



\B\/2wM eS all < 1. Defining = \^{z))e l ^ az -^ and 
\tp) = \(p{z))e lul ( az ~ tS> ~ %13 , we approximate — d 2 « oj 2 and 
d z ss iwa such that u 2 + -^jd 2 w 2uj(uj + i-d z ). Thus, 



/ 



i—d z 
a 



w p1 


iB v uja 



\ 2Af eH 



o 

iB x uja 
' 2AT.fr 



2aM c!{ 
iB x uj 
2aM c{{ 
™2 




We further simplify the mixing matrix by defining |7i) = 
\^{z))e- l P and \ip) = \(p(z)) e - l!3 with d z (3 = au^u). 
The magnetic field and frequency can be expressed in 
terms of their comoving values: Bo = a 2 B and ujq = aui. 
In the subsequent analysis we drop the subscript-zero 
notation for comoving magnetic field values. Thus we 
can write 



iad z — 




"eff 



2A/ off 2M cff 



— By 

2M otf 

3 2 

g m cff -a 

—5— —iO z a 

2ui 




= 0,(2) 



To solve this system of equations we must make various 
simplifying assumptions. We neglect spatial fluctuations 
in the electron density along the path length, and only 
consider a simple redshift scaling: n c (z) ~ ribo^eO" 3 , 
where X e is the ionization fraction along the path. The 
ionization fraction drops rapidly at recombination from 
X e ~ 1 to its final freeze-out value X e ~ 5 x 10~ 4 fljjj . 
and only increases again at the epoch of reionization. The 
conversion between photons and light scalar particles in a 
magnetic plasma scales inversely with the electron num- 
ber density in the plasma 12[ . Thus the conversion will 
be suppressed for high values of the ionization fraction. 
We assume the dominant contribution to photon-scalar 
mixing in a primordial magnetic field occurs between red- 
shift ~ 750 (by which time the ionization fraction has 
already dropped to 10 -3 ) and redshift ~ 20 (onset of 
reionization). We approximate the ionization fraction as 
being at a constant value of X e ~ 5 x 10~ 4 in this region, 
and neglect contributions from other sections of the path 
length. 

The chameleon mass for the generalized self-interaction 
potential of Eq. ([TJ is 



a 3 M 2a 4 M cS 



->+i 



8.4 x 10- 57 a- 9 / 2 ( ' GeV 2 



^10 9 GeV J 



where the second line assumes n = O(l), a magnetic field 
of less than 10 _9 G, a chameleon coupling strength of less 
than 10 _9 GeV _1 , and that the average baryonic density 
is 4% of the critical density. We compare this to the 
plasma frequency, 



4ira c 



4wa em X e p b0 



1.73 x 10 



-49-3 



m e m b 
GeV 2 



i p a 



over the redshift range ~ 20 to 750, where we have taken 
the average mass per baryon to be m& ~ 937 MeV [2(| . It 
is clear then that w pl 3> m 2 ^ along the path from recombi- 
nation to reionization. Our subsequent analysis requires 
the chameleon to be sufficiently light for the plasma fre- 
quency to dominate over the path length, and thus is not 
dependent on the specific form of the chameleon poten- 
tial. It applies to any scalar ALP satisfying this condi- 
tion. 

The comoving distance, z, travelled by the photons is 
related to the value of the scale factor at z by 



dz = da/a 2 H(a) . 



This leads to 



a(z) 



,1/2 



(0) + -H Q 



1/2. 



in the matter dominated era between recombination and 
reionization, where f2 m o = 0.26 [2l|. Over the total path 
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length from a redshift of 750 to 20, the comoving distance 
covered is L to t — 2.1 /i -1 Gpc. 

Changing the integration variable in Eq. (|2|) to 



-2 



such that = ad z , and noting pq/2luo 3> d z a, we find 

-By(tj) 








2M a 



2M ef£ 

Po 

2M e ff 2M eff 2w 




(3) 



where \x) = The ^-dependence of the mixing 

matrix is entirely due to fluctuations in the magnetic 
field. 



III. THE PRIMORDIAL MAGNETIC FIELD 

The generic model for the primordial magnetic field is 
of a stochastic field parameterized by a power-law power 
spectrum up to a cut-off scale ko, 

P(k)=A B k nB , k<k D , 

where Ab is some normalization constant. The spectral 
index must be greater than —3 to prevent infrared 
divergences in the integral over the power spectrum at 
long wavelengths. Assuming a statistically homogeneous 
and isotropic field, the power spectrum is defined by 

(B t (k)B*(k')) = (27r) 3 P lJ P(k)5(k - k') , 

where Py = Sij — kikj is the projector onto the transverse 
plane imposed by the divergence-free nature of the mag- 
netic field. We adopt the Fourier transform convention, 



Bi(k) - J d 3 xSi(x)e- 



•ik-x 



The two-point correlation function is then 

(Bi(x)fl 3 -(y)> = ^ J d 3 kP t3 P(ky k <*-^ 



Following the prescription in [5J, [22| , normalization of the 
magnetic field is achieved by convolving the field (in real 
space) with a Gaussian smoothing kernel of comoving 
radius Xb ■ Defining the Fourier transform of the Gaus- 
sian smoothing kernel to be /(k) = e~ x B k / 2 , the Fourier 
transform of the convolved field is 

Bi(k)Ua =Bi(k)-/(k). 

The mean-field amplitude of the smoothed field, B\, is 
given by 

2A B 1 /n B +3 



(2tt) 2 A 



iib+3 



Hence the normalization constant, 



A, 



(27T) 



n ; ■ "i Ij 2 



2fc™ B+3 r i 



"1+31 



where k\ = 2tt/\b- 

Constraints on the magnitude of the primordial mag- 
netic field from a comparison of Faraday rotation effects 
in the CMB with the WMAP 5-year data were given in 
@. They found that the upper limit on the mean-field 
amplitude of the magnetic field on a comoving length 
scale of As = IMpc was in the range 6 x 10 -8 to 
2xlO~ 6 G (95%CL) for a spectral indexes = -2.9 to -1. 
This range for the spectral index was based on the likely 
formation scenarios for the primordial magnetic field [2|1 
and current exclusion bounds on the spectral index [23j. 
More recent results analysing the WMAP 5-year data in 
combination with other CMB experiments @ constrained 
the mean-field amplitude -BiMpc < 2.98 x 10~ 9 G and the 
spectral index ns < —0.25 (95% CL). An analysis of the 
latest WMAP 7- year data Q derived upper bounds of 
SiMpc < 5.0 x 10~ 9 G and n B < -0.12 (95% CL). 

The magnetic field power spectrum will be damped at 
small length-scales. We assume a cut-off to the power 
spectrum at the Alfven wave damping scale [l5j]. This 
damping scale is approximated [22] as 



/ k L 



VMpc -1 



n B +5 



2.9xl0 4 ft 



B X V 2 / k x ^ nB+3 



10- 9 C; VMpc" 1 



where Hq = 100ft. km s _1 Mpc _1 . 

Using these definitions, we can calculate the correla- 
tions in the magnetic field along the line of sight: 

2?T 7T oo 

(Biix^Bjiyz)) = ^3 J J J k 2 sm9 k dkde k dcb k 



k 



where we have written the wave vector in spherical polar 
coordinates, k = (k,0k,<t>k)- Thus 



(B x (xz)B y (yz)) = (B y (xz)B x {yz)) = 0, 



and 



(B x (xz)B x (yz)} = (By(xz)B y (yz)) = R B (x - y), 
where 



Rb(z) 



(27T) 



k 2 P(k)dk 



sin(zfc) cos(zfc) sin(zfc) 
zk {zk) 2 (zk) 3 



Substituting the assumed form for the magnetic power 
spectrum, and changing variables, we find 



Rb{x -y) = 



B\K [(x - y)k D ] f2nk D 



Ilfl+3 



(4) 
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with 



and 



K[6] 



t nB+2 dt 



sin(0i) cos(6t) sin(0i) 



et 



{dtf 



(etf 



(5) 



IV. 



EVOLUTION OF THE STOKES 
PARAMETERS 



Following a similar procedure to that in jlq ] we subdi- 
vide the path length into many small domains of length 
L, each with a constant magnetic field strength and di- 
rection, and take correlations in the field across the dif- 
ferent domains to be described by the magnetic power 
spectrum. The inverse power-law form assumed for the 
primordial magnetic power spectrum imposes a rapidly 
decreasing amplitude to the field fluctuations and an even 
faster decreasing slope to the magnetic power spectrum 
with increasing k. At very large k the amplitude of the 
fluctuations are damped to zero by Alfven wave dissipa- 
tion. However slightly above the Alfven wave damping 
scale, k^ 1 , the amplitude of the field is non-zero yet the 
slope of the power spectrum will be approximately flat. 
In what follows we assume on scales slightly greater than 
kp 1 the magnetic field can be approximated as constant 
and take L ~ O (fc^ 1 ) • The exact value of L is only nec- 
essary in determining the number of magnetic domains 
that the photons traverse along the path. Some variation 
in L will not affect the order of magnitude of the final 
results. 



A. Evolution through a Single Domain 

For a single domain we define B to be the magnitude 
of the transverse component of the magnetic field, and 
<7 to be the angle it makes with the cc-axis. Within the 
domain these stay approximately constant. Returning to 
Eq. ([3]), we rotate the ("J x ,%) basis so as to reduce the 
problem to two-component mixing. Defining 

|%) = cos + sin a | %) , 
|7b) = -sma\%) +coscr|7y), 

we find 1 7 a ) is constant over the domain and 



idf 



( B/2M 
I B/2M -po/2wo 



\x) 



0. 



Diagonalisation of this two-component mixing leads to 

|7£ ew (L)) = |7°ew( O ^ e -i(A-A/co S 20) j 



|x ncw (L)) = |^ ncw (0))e _ '^ A+A / cos2e ^ 



where 



| -„cw ) 



cos 1 7f, ) + sin0|x) , 
-sin0|7 b ) +cos0|x) , 



tan 20 



2lj q B 
PqM cS " 



(6) 



A - -^KW-tfO))*-^. (7) 



Additionally defining 

A = sin 20 sin 
tan t/j = cos 20 tan 



A 



cos 20 
A 
cos 20 



the chameleon and photon polarization states in the ro- 
tated basis after passing through a single domain of 
length L, are 

\%{L)) = |7a(0)), 



|76(L)) = e 'Vl -A2|7 6(0))+ ? e^ A A|x(0)), 
\X(L)) = e - J Vl-A 2 |x(0))+*e- lA A|76(0)), 

where a = ip — A and (3 — ip + A. The probability of 
conversion between chameleons and photons in a single 
domain is P = A 2 . 

The intensity and polarization state of radiation is 
best described by its Stokes parameters: intensity, /, lin- 
ear polarization, Q and U, and circular polarization, V . 
These are defined in terms of the photon polarization 
states by 

I = (lx\lx) + (ly\ly) , 
Q = (lx\lx) ~ (jyhy) ) 

U + iV = 2( 7a; |7 y ). 

We additionally define chameleon 'Stokes parameters' to 
close the evolution equations, 

J + iK = 2e^( 7x |x), 
L + iM = 2e^( ly \ X ), 

i x - (xlx>- 

The evolution of the Stokes parameters through a sin- 
gle magnetic domain, in the rotated (%, %) basis, is then 



A 2 \ A 2 



Q 



u - 



+aVi - 


A 2 






)q 


-Ay/l- 


A 2 


e ia y/l- 


A 2 


+ie~ lP A 


(J- 



A 2 
2 



-7 * M X 



and 



L + iM 



J + iK -> e~ ll3 \/l-A 2 [J + iK 



+ie ta A(u + iV^j , 
e' 2 ^ (1 - A 2 ) (L + iM) 



+iA^l - A 2 (L-Q-2L 



I x is found by requiring that the total flux of photons 
and chameleons along the path length is conserved, I x = 
(hot ~ I 7 ) /a 2 . 



(SQ) (SU) 



(Qo) 



(Uo) 



-NR B (0) 



( 2oj 
\p M eS 

N-ln-l \ 

^ cos2A(n — r)i?B(z r — z n ) 



^2 I 2u; 



^^b(0)] 2 



N-ln-l \ 



n=0 r=0 



B. Evolution through Multiple Magnetic Domains 



The evolution equations above can be extended to 
many domains following a similar procedure to that in 

For CMB photons in the range 30-600 GHz, passing 
through a magnetic field of less than 10 _9 G, and assum- 
ing a photon-scalar coupling strength no greater than 
10 _9 GeV _1 , the mixing parameters defined in Eqs. ([5]) 
and ([7]) satisfy 6 <C 1 and | A | 1 in a domain of length 
L ~ hp 1 . There is a region towards the end of the path 
at larger a for which A ~ 1 , but the approximations that 
follow involving averaging terms in sin A and cos A along 
the path are still valid in this instance. Thus on average 
over the total path length, 



A 2 ~ 26 2 < 1 . 



(SV) 
(Vo) 



N—l n-1 

Y Y C0S 2 A ( n ~~ r ) i?B ( Zr 



n=0 r=0 



-2A 2 



PoM e ft 



N[R B (0)} 



where Rb(x — y) was defined in Eq (U). Over the range 
n B = -2.9 to -1, the integral K(9) (Eq. (J5J) can be well 
approximated by a window function: W(9) — 2/3(ris+3) 
for \6\ < 2, and zero otherwise. Taking the domain length 
L ~ fc^ 1 , the average modification to the CMB intensity 
is then 



(Sly) -NB 2 (2nk D /k x 
(Io) ~ 



,n B +3 



3(n B + 3)T(^^) \ Po M cS 



2uj 



, (8) 



This places us in the regime of weak-mixing which can 
be solved analytically if we neglect terms smaller than 
0(NA 2 ), where N — L tot /L is the number of magnetic 
domains along the path. Details of the mixing equations 
in multiple domains are presented in the Appendix. Here 
we quote the average, over many lines of sight, of the 
modification to the Stokes parameters. We assume there 
is no initial flux of chameleons and define z n ~ nL to be 
the location of the n th domain. Then, 

N-ln-l \ 
+ Y Y C0S 2 A ( n ~ r ) i?B ( 

Z'r Z n J I , 

n =Q r=0 / 



where we have averaged terms in cos A and sin A, since 
|A| 3> 1. Of the average modifications to the Stokes pa- 
rameters, that of the CMB intensity will be most sig- 
nificant since the CMB is only weakly polarized with 
(V 2 ) 1 / 2 «(Q 2 }^ 2 ,(US) 1/2 «Io. 

In addition to the average modification to the Stokes 
parameters as a result of photon-scalar mixing, there 
will also be a change to the correlations between dif- 
ferent Stokes parameters. A full derivation of the cross- 
correlations along different lines of sight would require us 
to expand the primordial field power spectrum in spher- 
ical harmonics, a task we reserve for a future paper [24j |. 
However we can estimate the magnitude of the effect by 
looking at, say, the (UV) cross-correlation along the line 
of sight. The dominant contribution, neglecting terms 



8 



proportional to (QqUq}, (IqVo) and smaller, is 



N-l 



\ n : m—0 

N-l B-lm-1 

S ' cos 2 A (ri-r) sin 2A(m- s)R B (z n - z m )R B (z r - z 3 ) 



n,m— r— s=0 
JV-1 n-l 



+ A ^ y^fl B (^ n -Zm)i?.B(z r -Zm)cos2A(n-r) (JqQq 



n,m— r— 



-2AA 



5^(27Tfcp/fc A )"B+ 3 

(n B +3)r(^|±^) 



(ifc)VoQo). 



From this we see that a significant (E/V) correlation can 
arise from the much larger (IQ) correlation. Although 
the above analysis only applies to the average correla- 
tions along many lines of sight, for which (IqQq) ~ 
in the CMB, we expect a similar magnitude to the effect 
when a full derivation for the cross-correlations along dif- 
ferent lines of sight is made. It opens up the exciting 
possibility of a significant chameleon contribution to the 
(EB) power spectrum of the CMB. 



V. BOUNDS ON THE CMB INTENSITY FROM 
FIRAS 

Returning to Eq. (|5J|, we have that the average modi- 
fication to the CMB intensity from scalar-photon mixing 
in a primordial magnetic field is 



(Sly) _ ~ 
— l 7 _ 

1 y0 

where we have defined 



\p M e g 



B, 



4 Bl(\ B k D ) n » +3 
3 (n B + 3)T(^) 



This expresses the primordial magnetic field in terms of 
its mean- field amplitude, B\, at a length-scale of A^. 
Most constraints on the magnitude of the magnetic field 
are expressed in terms of these parameters. The value of 
us dictates the slope of the magnetic power spectrum, 
and we consider the range —2.9 to —1. 

The most precise measurements to date of the CMB 
monopole near the peak in its spectrum come from 
the far infrared absolute spectrophotometer (FIRAS) on 
board the cosmic background explorer (COBE) satellite 
[25l [26| . These measurements fit exceptionally closely 
to the spectrum of a black-body at a temperature of 
2.725K ± lmK. In Table 4 of [H, the residuals of the 
CMB monopole that remain after subtracting foreground 
effects are listed. These residuals, added to the spec- 
trum of a perfect black-body radiating at 2.725 K, can 
be compared to the predictions of photon-scalar mixing 
to constrain the parameter space. The FIRAS instru- 
ment has 43 significant frequency channels in the range 



60-630 GHz after calibration. The spectral intensity of a 
black-body at temperature To is 



lo 0, To) 



4iris 3 



exp {2-kv/Tq - 1) 



In the presence of photon-scalar mixing the spectrum is 
modified to 

I(u,T ,X) = (l-P 7 _^)Jo(i/,T ) , (9) 
where A is defined such that 

^^.TAXIO-G^) 2 , 



A (B e g, M e g) — 



B eS V ( M eff V 



io- 9 Gy V 109GeV / 

We maximize the likelihood, L, defined by 



-21ogL = ]T 



43 /job 
i=l 



fT° fcS (^) 



over the parameter space (To, A). Here I obs is the ob- 
served CMB monopole measured by FIRAS at 43 differ- 
ent frequencies, 1^, with standard errors a obs . Tq, A) 
is the spectral intensity predicted for chameleon-photon 
mixing given in Eq We treat To as a free parameter 
in our analysis rather than fixing it at Tf IR as = 2.725 K. 
Confidence limits are estimated by assuming 



X 2 = - log 



L(T Q ,X) 



follows a Xi distribution, where ^To,A^ are the best fit 

values of To and A in the parameter space, found by the 
maximum likelihood procedure. Minimizing \ 2 with re- 
spect to A gives T = 2.725 ± 0.016 K, identical to the 
value obtained when no chameleon mixing is considered. 
Minimizing with respect to To gives the following 95% 
confidence limit on A, 

log 10 A < -0.01 (95%). 

Fig. [j] shows the 68%, 95% and 99.9% confidence limits 
on the (To, A) parameter space. 

The limit on A corresponds to a degenerate constraint 
on the magnitude of the primordial magnetic field and 
the chameleon-photon coupling strength. The magnetic 
field strength is expressed in terms of the mean-field am- 
plitude at a comoving length-scale of 1 Mpc, B\. In Fig. [5] 
we plot the exclusion bounds in the (Bx^M^g) parame- 
ter space resulting from the 95% confidence limit on A. 
Different lines correspond to different values of the mag- 
netic spectral index, ns, in the range -2.9 to -1. The 
open squares mark the 95% upper limit on the allowed 
values for B\ (for each ng) found in Q from compari- 
son of Faraday rotation effects in the CMB with WMAP 
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FIG. 1: Confidence limits on the black-body temperature, To, 
and the strength of photon-scalar mixing, parameterized by A, 
coming from precision measurements of the CMB monopole, 
where Tfiras = 2.725 K. 
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FIG. 2: Exclusion bounds on the primordial magnetic field 
and photon-scalar coupling strength, g c R = 1/M e g, from pre- 
cision measurements of the CMB monopole, at the 95% con- 
fidence level. B\ is the mean-field amplitude of the magnetic 
field at a comoving length-scale of 1 Mpc. Open squares cor- 
respond to the upper limit on the primordial magnetic field 
found in [5J, while the vertical line corresponds to the upper 
limit found in @, 0] . The region to the bottom right of the 
plot is excluded for different values of the magnetic spectral 
index nj. 



5-year data. The vertical line corresponds to the more 
recent constraints on B\ given in Q and Q- The region 
to the bottom right of the plot for each rig value is ex- 
cluded at the 95% confidence level. For example if we 
take a magnetic field strength at the upper limit allowed 
by Faraday rotation effects in the CMB, B\ ~ 5 x 10 _9 G, 
the corresponding bound on the photon-scalar coupling 
strength, g cS = l/M off , is 

M eff > (0.14-1.09) x 10 13 GeV, 

depending on the slope of the magnetic power spectrum. 



VI. CONCLUSIONS 

The existence of a large-scale cosmological field in the 
early Universe is so far unconfirmed. It would need to 
be in the region of O (10~ 10 ~ 10~ 9 G) if it is to explain 
the formation of the O(fiG) magnetic fields in galaxies 
and galaxy clusters, through adiabatic collapse. If this 
primordial magnetic field exists it would induce mixing 
between CMB photons and axion-likc particles (ALPs) as 
they propagate from the last-scattering surface to Earth. 

In this paper we have studied the case of non-resonant 
mixing between scalar- ALPs and photons in a primordial 
magnetic field, with specific reference to the chameleon 
scalar field model. 

The chameleon model is a promising candidate for the 
dark energy scalar field since it can have a gravitational 
strength (or stronger) coupling to normal matter while 
at the same time evading fifth-force constraints. To date, 
the strongest bounds on the chameleon coupling strength 
come from chameleon-photon mixing in local astrophys- 
ical environments, such as starlight propagating through 
the galactic magnetic field: g e g < 9.1 x 10~ 10 [1(|. Should 
there be a detection of a primordial magnetic field of or- 
der > 10~ 10 G, our results would place far greater con- 
straints on the coupling strength. 

We have considered a stochastic primordial magnetic 
field described by a power-law power spectrum, P(k) oc 
k nB , up to some cut-off damping scale fcjj . Photon- 
scalar mixing in this field was solved by dividing the path 
length into multiple magnetic domains in which the field 
is approximated as constant. The length of the domains 
was taken to be of a comparable size to fc^ 1 since the 
magnetic power spectrum will be approximately flat, but 
non-zero, just above the damping scale. Correlations be- 
tween the magnetic field strength and direction in each 
domain are determined by the magnetic power spectrum. 
In addition, we approximated the ionization fraction as 
being constant in the region from redshift ~ 750 to 20, 
and held at its post-recombination freeze-out value of 
X e ~ 5 x 10~ 4 . The dominant contribution to photon- 
scalar mixing in the CMB will occur in this region of 
low electron density, and we neglected contributions from 
other elements of the path length. A more sophisticated 
approach to modelling the electron density from recom- 
bination to the present day may result in small changes 
to the predictions, but would require a numerical rather 
than analytical approach to the mixing equations. 

We have compared our predictions of the average modi- 
fication of the CMB intensity over the whole sky, to preci- 
sion measurements of the CMB monopole by the FIRAS 
instrument on board the COBE satellite (2f| [2(|. This 
constrains the probability of photon-scalar mixing over 
the path length, to be 

-P 7 ^ < 0.026 (95% CL) 

at 100 GHz. The corresponding bounds on the magni- 
tude of the magnetic field and photon-scalar coupling 
strength are plotted in Fig. [2] for different values of the 
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magnetic spectral index. Until a detection of the primor- 
dial magnetic field is made, to break the degeneneracy of 
this constraint, we cannot place limits on the chameleon- 
photon coupling strength, g c g. For the largest magnetic 
field allowed by the constraints in [y, |7j we would find 
the strongest possible constraint, 

g cS < (0.92 - 7.14) x 10" 13 GeV _1 , 

depending on the slope of the magnetic field power spec- 
trum. 

The results in this paper apply to any scalar ALP with 
a mass less than ~ 10 _14 eV, since we require the mass 
to be much smaller than the plasma frequency along the 
path. These nicely complement the bounds derived in 
[1J| for resonant conversion between photons and ALPs 
in a primordial magnetic field, g e s{B 2 ) 1 ^ 2 < 10~ 13 ~ 
10~ n GeV -1 nG, which apply to ALP masses in the 
range 10" 14 eV to 10~ 4 eV. 

In addition to the average modification to the CMB 
intensity, the formalism developed in this paper can 
straightforwardly be extended to calculate the change to 
correlations between the CMB Stokes parameters along 
different lines of sight. In section IIVB1 an example 
was given of how a correlation between the U and V 
polarization modes can arise from photon-scalar mixing 
in the primordial magnetic field, given a non-zero (IQ) 
correlation. If this effect is present in the CMB cross- 
correlations, it would lead to a significant chameleon sig- 
nature in the (EB) power spectrum. A full analysis of 
this effect is kept for a separate publication (24|. 



relations for the Stokes parameters: 



(51 



J J n+1 



SQ 



5V n 



A 2 A 2 

(M 7 ) B - -^J - -y (cos 2ct„Qo + mn2a n U ) 

+A n cos a n [L n sin 2 A — M n cos 2 A] 
+A n sin a n [K n cos 2A — J n sin 2 A] , 



n+1 



SQ r 



A 2 



-Qo 



A, 



• cos 2<7„/o + a n sin 2a n V n 



— Y sin 2a n (sin 2a n Qo — cos 2a n Uo) 

—A n cos a n [L n sin 2A — M n cos 2A] 
+A n sin cr„ [K n cos 2A — J n sin 2 A] , 



SU n 



A', 



-U 



■ sin 2<j n 7o — en cos 2a n V n 



— Y cos2er„ (cos2ct„[/o — sin2cr„Qo) 
— A n sin a n [L n sin 2A — M n cos 2A] 
— A n cos a n [K n cos 2 A — J n sin 2 A] , 

A 2 a 2 
SV n - -fV - -fV 

+a n (cos2cr„f7„ - sin2cr n Q„) 
+A n sin cr„ [L n cos 2 A + M n sin 2 A] 
+A n cos a n [K n sin 2 A + J„ cos 2 A] , 
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Appendix A: Photon-Scalar Mixing in Multiple 
Magnetic Domains 

In section lTV Al we found the evolution equations of the 
Stokes parameters after passing through a single mag- 
netic domain of length L, within which the magnetic 
field strength and direction are approximated as con- 
stant. This can be extended to N multiple domains, 
under the assumption of weak-mixing for which we re- 
quire NA 2 < 1 and N > 1. CMB photons in the 
range 30-600 GHz, propagating through a weak primor- 
dial magnetic field, satisfy 6 <C 1 and |A| ^> 1. In this 
limit ij) ~ A, ~ 2A and a - O(A) < 1. We de- 
fine (5/n+i to be the difference in intensity after passing 
through the n th domain compared to its initial value, Jo, 
and similarly for the other parameters. Neglecting terms 
smaller than 0(A 2 ), we obtain the following recurrence 



J„+i + iK n+ i ~ e~ 2lA ( J n + iK n ) - A n cos a n (V - iU ) 

-iAnSman (I + Q ) , 
L„+i + iM n+1 ~ e~ 2lA (L n + iM n ) - A n sin a n (V + iU ) 

+iA n cos a n (I - Qo) , 



where we have assumed there is no initial chameleon 
flux which necessarily sets \x) ~ Note that A 

depends on the magnitude of the transverse component 
of the magnetic field, which fluctuates over the different 
domains, hence the subscript n, whereas A is constant 
along the path save for a weak dependence on the scale 
factor. However, since |A| 3> 1, this will be averaged out 
across the path length. Solving this system of equations 
for propagation through N domains, we find the modified 
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Stokes parameters: 



. N-l n-l 



(i - \v N ~ r N ) I + {\k c n + &*) Q 



Q ~ (1 - |P W - A 2 [iyyf + ^] - 0£r) Q 



A 



2^ + Ma? 



N 



U 



U 



+ + *o + KA + e#) Vo , 

(i - \v N - a 2 + n%] - r N ) u 



2 W + A% 



^" ) Go 



V" 



+ QkJt + j o - («atA + g c +) V , 

(i - \v N - a 2 [iQjy + - - r N ) V 

- (k s n A + g s ^)Q + (k c n A + g c +) U - g s N I , 



where we have defined 
^N-l 



= Y,n=0 A n > = E!Lo' A n ) 



and 



= J2n=o^2r^o A n A r cos (2A(n- r)) cos (a r ±a n ) , 

$ S N = ^n=0^2r=0 A n A r COS (2A(n-r)) Sill (a r ±a n ) > 

Q C N = Efjo 1 E"=0 j4 " i4 <- Sin ( 2A (™ ~ r )) C0S ± ^-i) ' 

Q S n = J2n=a S"=o A n A r sin (2A(n - r)) sin (a r ± a n ) , 



and 



= T,n=oT^=0 A l A r COs2 ^rCOS2a n , 

= Y,n=oY,r=0 A l A l sin2 < T rCOs2(J n , 

= E^oE^g 1 A 2 A 2 cos2CT r sin2CT„, 

t^N = 12n=oT,rZa A 2 A 2 sin2a r sin2cr„, 



and 



K Af — En=0 A n c OS2CT„ , K% — Yln=0 Sm 2<7„ , 

»?" = E^o 2a„ , ^ = En=o 4 cos 2 2a„ , 

Vn = Y,n=0 A n sin 2cr ™ COS 2cr " ■ 

We define z n to be the location of the n th domain, then 
B x {z n ) = B n cosa n and B y (z n ) = B n sin <r n by defini- 
tion. Remembering A n = LOoB n /poM e s, we find that the 
average over many lines of sight of the above quantities 
are 

<^> - 4 ^fe) 4 [^(o)] 2 , 



(Mat s > 



p M e j 



M ^2 J2cos2A(n- r)R B {z r -z n ). 



n=0 r=0 
, AT-1 n-l 



n=0 r=0 

4 AT-1 n-l 

,i=0 r=Q 

(^)= iV (iafe) 4 [ i? B(0)] 2 , 



and all other quantities average to zero. We have as- 
sumed fluctuations in the magnetic field are approxi- 
mately Gaussian so that the four-point correlations can 
be expressed in terms of the two-point correlation func- 
tion, R B (x — y), defined in Eq. Qj. 

Taking the average of the modified Stokes parameters, 
over the whole sky, we find 



(«7 
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2ojq 
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JV-ln-l \ 
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(SQ) (SU) 
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(Uo) 



\poM eS J \2 BV ' 



N-l n-l 



2_\. cos 2 A(rt — r)i?B (Zr — z n ) 



(SV) 
(Vo) 



2uj 

N-l n-l 
n=0 r=0 



/ 2w 



\n[r b (o)} 2 



\p M eS 



-NR B (0) 



N-ln-l \ 

^2 E C0S 2 M n - r ) R B (z r - Zn) 
n =0 r =0 ) 



These results are quoted in section HVBl 
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